ENERGY MINIMIZATION, PERIODIC SETS AND SPHERICAL DESIGNS 

RENAUD COULANGEON AND ACHILL SCHURMANN 

Abstract. We study energy minimization for pair potentials among periodic sets in Eu- 
clidean spaces. We derive some sufficient conditions under which a point lattice locally 
minimizes the energy associated to a large class of potential functions. This allows in 
^^ ' particular to prove a local version of Cohn and Kumar's conjecture that A2, D4, Eg and the 

fSJ ! Leech lattice are globally universally optimal, regarding energy minimization, and among 

f^ • periodic sets of fixed point density. 
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^ ■ 1. Introduction 

. 5r , The study of point configurations that minimize energy given by some pair potential 

occurs in diverse contexts, such as crystallography, electrostatics or computer graphics. 
There exist numerous numerical approaches to find locally optimal or stable config- 
urations. However, a mathematical rigorous treatment proving optimality of a point 
configuration is quite difficult. 

Already in 1897, J.J. Thomson, the inventor (discoverer) of the electron, came to the 
conclusion that "the equations which determine the stability of such a collection of par- 
ticles increase so rapidly in complexity with the number of particles that a general math- 
ematical investigation is scarcely possible". In some special situations more can be said 
though. 



2000 Mathematics Subject Classification. 82B, 52C, iiH. 

Key words and phrases, energy minimization, universal optimality, periodic sets. 



1.1. Energy minimizing spherical codes. One important case that has been studied ex- 
tensively by physicists are point configurations (charged particles) on the surface of a 
sphere. We call such sets spherical codes in what follows. We may consider the unit 
sphere S^~^ C R^. Given a real- valued nonnegative function / : (0,4] — ^ R, we ask in 
this situation to minimize the /-potential energy 

w Eif,c)= E /(ll^-yf) 

among point sets C on S'^^^ of fixed cardinality |C|. 

One may think of C as a set of electrically charged particles in R^^, and of the potential 
/(r) = r^^^, in which case definition ^ coincides with the classical notion of potential 
energy in physics. Depending on the function /, the solutions may look quite differ- 
ent. However, for many "reasonable functions" /, we may get the same solution. For 
example, the vertices of a regular simplex (tetrahedron in R'^) will be optimal for any 
continuous and decreasing function /. 

In | |CKo7[ , Cohn and Kumar introduced the notion of a universally optimal configuration 



of points. On the unit sphere, they are point configurations C that minimize E(f,C) for 
all completely monotonic functions, that is, for all real- valued C°° functions on the interval 
I = (0,4], such that (-1)*'/^''^^) > ^r all x in 7 and all k>0. The class of completely 
monotonic functions contains in a sense all the "reasonable functions" in the context of 
energy minimization. It in particular contains all the inverse power laws f{r) = r~^ with 
s > 0, which are often studied in physics. 

It turns out that there are several fascinating examples of universally optimal spherical 
codes. Indeed, Cohn and Kumar were able to derive a sufficient criterion for universal 
optimality, via so-called linear programming bounds and they showed that all spherical 
configurations with at most m mutual distances, that form a spherical (2m — 1) -design 
(see Definition [472]), are universally optimal. This shows by example that there are excep- 
tional structures (and infinite families of them) for which a rigorous mathematical proof 
of a very broad energy minimization property is possible. 



1.2. Energy minimizing periodic point sets. Cohn's and Kumar's considerations for 
spherical point sets can be extended in several ways. Note first that definition ((Tj makes 
sense for any finite set in R'^. One further natural extension is to consider infinite sets 
C in Euclidean spaces, possibly unbounded. This however yields difficulties in defining 
potential energy properly, because of possible subtle convergence problems. For periodic 
sets such problems can be avoided. We say a discrete set in R^ is periodic, if it is a disjoint 
union of finitely many translates of a given full-rank lattice L C R'^ (a discrete subgroup 
of R'*). In particular, a full-rank-lattice L C R'* itself is a periodic set. In general we can 
write. 



A=\Jit, + L) 
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where ti, . . . , tm are some vectors in R . For a potential function /, the /-energy of A is 
defined as 

-1 m 

(2) £(/'A) = -E E f{\\i^-A?)■ 

It can be shown that, when finite, the right-hand side of O is equal to 

(3) lim r^ E m^-yf) 

I ^^1 x,yeA[i,xty 

where Ar = {x G A, ||x|| < R} (| CKo7[ Lemma 9.1.]). This clarifies the link with the 
definition of energy for finite configurations of points. 

As in the case of spherical codes, one may ask if there exist universally optimal periodic 
sets, that is, periodic sets that minimize the energy E(/, A) for all completely monotonic 
functions /. At this point, no such universally optimal periodic set is known. However, 
exceptional structures as the hexagonal lattice, the root lattice Eg, and the 24-dimensional 
Leech lattice are conjectured to be examples (see | CKo7| ). Recent experiments show (see 
[CKSo9 1) that also the root lattice D4 and (somewhat surprisingly) the periodic non- 
lattice set Dg could be universally optimal. 

As a first attempt to prove universal optimality for any of the examples above, it is 

natural to ask whether universal optimality holds at least locally. Before we go further, 

we recall a few known results for a similar question in the noticeably simpler context 

1 
of lattices. Indeed, when A = L is a lattice, and /(r) = — for some s > j , then the 

1 
corresponding energy E{—,L) coincides with the Epstein zeta function of L 



Similarly, if f{r) = e~^^ , the corresponding energy is 

c 
n 



£(e— ,L)= E e-'\\-\\=e^{i-)-l, 
o^xeL 



where 61 is the usual theta series of L. Questions of optimality for lattices with re- 
spect to their zeta (resp. theta) function have recently been investigated by Sarnak and 
Strombergsson in IISS06I , and by the first author in IIC0U06I , in connection with the the- 
ory of spherical designs. In particular, one has the following sufficient condition for local 
optimality among lattices: 

Theorem 1.1 ( |Couo6| ). 

Lattices for which all shells are ^-designs achieve a local mimimum (among lattices) of the map 
L I— > E{e~'^\ E) for big enough c. 

Cohn and Kumar observed in ]CKo7[ that it is enough to deal with potentials of the 
type r 1— > e~'^^, c > 0, to recover all completely monotonic potentials. So, in view of 
the above theorem, universal optimality among lattices essentially reduces to a prop- 
erty of theta series. Nevertheless, to actually infer universal local optimality (among 
lattices) from Theorem [iTTl one has to be able in addition to remove the restriction to 
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"big enough c" , and get the result "for any c > 0" instead, which is highly non trivial 
in general. It turns out to be possible however in the case of A2, D4, Eg and the Leech 
lattice, thanks to Sarnak and Strombergsson's result in |,SSo6,l (see IISSo6[ Theorem 1]). 

As in the case of spherical point sets, where the kissing number problem can be seen as 
a limiting case, the sphere packing problem (asking for the maximum possible minimum 
distance of points at a fixed point density) is a limiting case of energy minimization 
of Euclidean point sets. It can be shown that the density of periodic packings come 
arbitrarily close to the optimal density of a sphere packing in a given dimension d. 
Whereas the local optima for the density of lattice packings are well understood through 
Voronoi's characterization in terms of perfection and eutaxy, the situation for periodic 
packings is comparatively more difficult. It was shown in HSchiol that if a lattice L is 
perfect and strongly eutactic {i.e. the minimal vectors form a 2-design), then L achieves a 
local maximum for density not only among lattice sphere packings but also among all 
periodic sphere packings HSchioi Theorem 10]. By a theorem of Venkov, the condition 
that L be perfect and strongly eutactic is satisfied in particular when the set of minimal 
vectors forms a 4-design. Lattices satisfying this property are sometimes called strongly 
perfect in the literature. To summarize, one has 

Theorem 1.2 ( HSchiol ). 

Lattices for which the set of minimal vectors forms a ^-design achieve a local optimum for the 
sphere packing problem among all periodic sets. 

In this respect, strongly perfect lattices are somehow extremely rigid : there is no 
possibility to improve locally their density within the set of periodic sets. 

The aim of this paper is to combine the ideas of IIC0U06I and HSchiol to prove essen- 
tially that lattices satisfying the conditions of Theorem lTTTI are locally universally optimal 
not only among lattices, but indeed among all periodic sets (a precise formulation is 
given in Section [4J Theorem ^^ . Again, this means that the 4-design property yields 
a strong rigidity. Our main result (Theorem [44]) can be in particular applied to the lat- 
tices A2, D4, Eg and to the Leech lattice (see Theorem \/[.6) : this generalizes the result of 



Sarnak and Strombergsson in IISS06I , and proves a local version of Cohn and Kumar's 
conjecture [ |CKo7[ Conjecture 9.4]. 

2. Preliminaries. 

2.1. A space of parameters. The study of local variations of energy first requires a suit- 
able parametrization of the space of periodic sets. From now on, unless otherwise stated, 
the word "lattice" will stand for "full-rank Euclidean lattice", i.e. for a discrete sub- 
group of maximal rank in the Euclidean space M.'^, equipped with its standard norm 

1 



\Ili=i ^j ) ■ Following | Scho9[ we say that A C R'^ is an m-periodic set if there 



exists a lattice L C IR^ and vectors t\, . . . ,tyn in IR'^ such that 

m 

(4) A = LJ {U + L), 

z=l 
4 



the disjoint union of m translates of L (in other words, we assume that the sets ti + L are 
pairwise disjoint, i.e. ti — tj i L for i ^ ;). We denote by Cm the set of m-periodic sets in 



R'^. 



We define the point density of A = \JIL\ U + i- as 



m 



This accounts for the number of points per unit volume and is of course independent 
of the representation of A as a union of translates of a lattice (we use the terminology 
point density rather than simply density to avoid any confusion with the density of the 
associated sphere packings). 

Since most of the quantities we will be considering {e.g. energy, packing-density) are 
invariant under orthogonal transformations and translations, we may identify two m- 
periodic sets which are isometric. In particular, the m-tuple {ti,...,tm) can be defined 
up to translation of its components by a common vector. In what follows, we adopt the 
notation 'Rf'^ to refer to the set of m-tuples u = (ui, . . . ,Um) of vectors in R'* subject to 
the condition 

(6) Ui - Uj i Z.^ for i i^ j 

and we denote by R'"'^/T the same set up to translation. For any u = (wi, . . . , Um) G IRf^, 
we define a standard periodic set 

m 

(7) nu = [J(u, + z^). 

Then, any m-periodic set may be written as AOy for some A G GL^(R) and u G R'^'^. 
The matrix A in the above expression is determined, up to left multiplication by 0{d), 
by the positive definite quadratic form Q = A'^ A. Note that we use column vectors, and 
with these settings have 

||Ax|p = Q[x\ ■.= x^Qx. 

Using the notation S for the set oi d x d real symmetric matrices and Sz^q for the cone 
of positive definite ones, we thus get a parametrization of 0(d)\£m/T by <S^q x R"''^/T: 
to (Q, u) G *S^Q X R|f^/T one associates the m-periodic set A = AOu, where A is a 
square root of Q. In keeping with |Scho9| l, the elements of S^q := S'^q x R"''* are called 



m-periodic forms . 

Finally, energy comparison between two m-periodic sets makes sense only if they are 
assumed to have the same point density (otherwise, by shrinking/expanding a given 
periodic set with a scaling factor, one can achieve arbitrarily small/large energy). One 
can for instance restrict to m-periodic sets of point density m, which amounts, in the 
above parametrization by periodic forms, to consider the space V^q := V^q x M.I^'^ /T, 
where V^q stands for the set of positive definite quadratic forms of determinant 1. 

In accordance with formula (Q or ([3]), computation of energy involves evaluating po- 
tential functions over the set of nonzero elements in 

A-A:={x-y : x,y e A} . 
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One difficulty is that a given element in A — A generally admits several representations 
as a difference of two elements in A. The situation is somewhat simpler when A is a 
lattice, as shown by the following lemma. 

m 

Lemma 2.1. Let A = \_^ tj + L be an m-periodic set in R^. For x G A, set 

i=\ 

Ax = {y - X : y G A} . 

Then the following assertions are equivalent 

(1) A is a lattice. 

(2) A - A = A. 

(3) Ax = A /or all x G A. 

(4) For any k in {1, . . . , m}, there is a uniquely defined permutation a^^ 0/ {1, ... , m} such 
that 

Vz G {1, . . . , m} fp.^(;) = ti - ffc mod L. 

Proof. Lattices are characterized as discrete additive subgroups of IR and the equivalence 
of dl]), ^ and (J3]) is derived from that. As for ^ =^ ^, we have that for fixed k, the 
difference ti — tj^ lies in Afj. = A = UjLi ^j + ^/ so there exists a uniquely determined 
index (7jt(z) such that f, — f ^ G fcr^(;) + L. Moreover, cr]^{i) = dkij) i^ ^rid only if ti — tj^ = 
tj — tj^ mod L, which means that ti — tj G L, whence z = /, so a^^ is a bijection. Finally, 
property (J4J) cleary implies that any pairwise differences of elements in A are in A, which 
shows that ([4]) ^ ©. n 

2.2. Potentials. As regards the potential functions / to be used, the following assump- 
tions will be made throughout: 

Assumption 1. / is a completely monotonic function on (0, 00), i.e. real-valued C°° functions 
on (0,00) such that (—l)^f^^\x) > for all x in (0, 00) and all k > (in particular, f is 
nonnegative). 

Such functions will be referred to as CM functions in what follows. 

It will also be useful, although not necessary, to assume that formula (tlj converges, 
which is ensured by the following 

Assumption 2. There exists e > such that f{x) = 0(x)~2~^ as x tends to infinity. 

It would be possible, using Bernstein's theorem, to restrict to potentials of the form 
/c(r) = e~'^'' with c > 0. Indeed, any CM function may be written as 



/■OO 

f(x) = / e-'^'daic) 
Jo 



(Stieljes Integral) for some weakly increasing function a. (see |Wid4i [[Theorem 12b, p. 
161]). 

A case that we consider separately on its own first, is that of inverse power laws 
Ps{r) = r~^ for some s > 0. These do not encompass the whole class of CM functions, 
but they are easier to deal with. 
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3- Local Study of Potential Energy. 

3.1. Local expression for the energy. The /-energy of an m-periodic set A = AOu 
depends only on the associated periodic form (Q, u), namely one has 

-1 m 

(8) E(/,A) = E(/,(Q,u)) = -E E /(QK-x]). 

We want to expand the /-energy in a neighbourhood of a given m-periodic set 

Ao = AqO^o = U^li^o + Lo, 

where we set Lq = AqZJ^ and t*^ = Aqvl^ {i.e. t^ = Au'f for 1 < i < m). We also assume 
that Aq has point density p^(Ao) = m and we let Xq = (Qo, u*^) be the corresponding 
periodic form, with Qq = AqAq. 

The manifold V^q = 7^>q x R|^'^/T is locally homeomorphic in a neighborhood of 
^0 = (Qo/ii'^) to its tangent space at Xq which is identified with Tqq x R"^'^/T where 

TQ, = {KeS' : Tr(Qo-iK) = 0} . 

The isomorphism is obtained via the matrix exponential through the map {K,u) 1— )• 
( Qo exp(Q^ K),u^ + u]. Note that the tangent space Tqq x R'"^/T at Xq comes equipped 
with its standard SL(i(R) -invariant scalar product 

m 
i=l 

which defines the Riemannian structure of V^^. To study the local variations of the 
/-energy around Xq, it is enough to consider the /-energy of {Qoexp(QQ^H),u^ + u), 
for small enough H G Tq^ and u G R'"'^/T. It equals 



-1 m 

(10) -E E /(Qoexp(Qo-iK) 



1=1 XEil n 



U^ + Ui — X 



Each term u^ + u, — x in the internal sum may be written as u'^ — u^ + w, — Uj + v for 
some ; G {1, . . . , m} and some i; G Z'^. Note that the condition u^ — u^ + ui — Uj + v i^ 
will be satisfied if and only ii u'^ — u^ + v itself is non-zero, provided that the w; are close 
enough to (this will be the case for instance if the Uj's satisfy ||uf || < ^, where po := 
mino^xen o-n ll-*^!!)- Consequently, assuming that u lies in a suitable neighbourhood of 
0, we can rewrite ((To|| as 

(11) ^ E E f{QoexpiQ^'K)[w + u,-Uj]) 

\<i,i<m Qtweu^-u°:+'Z<^ 



In order to get simpler expressions in the calculations to come it is more convenient 
to change coordinates, that is we rewrite the above expression as 

(12) £;(H,t):=l E E f{exp{H)[w + t,-tj]) 

l<i,i<m o^ii;6f°-t°+Lo 

where t = Aqu (resp. 1° = AqvP), and H = (Aq^) i^A'^ is now in 7Id = {H G 5^^ : Tr(H) = O}. 
To see that ((Ti|) and ((12)) are the same, we use the identity 

Qoexp(Q-iK)[x] = x'Qoexp(Q-iK)x = x^A[,exp Ua'^'^' KAq''\ Aqx = exp(H)[Aox] 

valid for any x G R^. That H = ( Aq^ j KAq^ is symmetric is clear, and Tr H = follows 
from the simple observation that 

Tr ((Aq-iJ^KA-i) = Tr (a"! {a-^)' K^ = T^Qq-^K). 

Note that the scalar product ^ on 7id/ which we denote simply by ( , ) in what follows, 
takes the form 

m 

(13) ((K,u),(L,v)) =Tr(KL) + Et/[z;, 

z=l 

Note that the definition of Ef depends on a given representation of Aq as a periodic set, 
that is, it depends on Aq and u*^. Note also that E(/, A) = £^(0, 0) with this setting. 

The two main ingredients to obtain further simplifications in the above formula are 
the following : 

(1) use the additive structure of Aq (if any). 

(2) use translation invariance of the energy. 

These conditions are met in particular when Aq is a lattice, in which case we obtain 
the following crucial lemma : 



m 
^0 



Lemma 3.1. Suppose that Aq = |_J (^^ + Lq) is a lattice. Then 

i=l 

-1 

(14) Ef{H,t) = -^ E E f (exp{H) \w + 1, - 1,^(^,^ 

where a> is the bijection of Lemma 12m 



Proof. Since Aq — Aq = Aq, any coset t^ — t^ + £0 iri the internal sum ([12)) can be written 
as t^ + Lq for a uniquely defined k. More precisely, using Lemma i2.iC [2J), we obtain 

-1 m m 

(15) E/(w,t) = -E E Lf[^MH)[^^ + u-t,,ir)\) 



where O/^ is the permutation defined by the condition that fi , .. = t'^ — t'^ mod Lq for all 
z G {1, . . . , m} (see Lemma im)). Note that the tj are replaced by t^ n-j and that the change 
from index / to k causes a reordering of terms. 

Because of the translation invariance of the energy, the energy is not modified if all the 
components of t*^ are translated by a common vector a. G R'^. In particular, we can choose 
a = —fi for some ; G {1, . . . , m}. Applying this to < ^T^ , we get for any ; G {1, . . . , m}, the 
equation 



Ef{H,t) 



-E E Ef[^MH)[^ + tr-t,,ir 



k=l Qi=ivetl-t°+Lo '=1 



m 



E E Ef(eMH)\w+u-t^^^, 



fc=i o^^,e-tl^^.^+Lo '=1 



Adding up the (Ej) ? for / = 1, . . . ,m and then averaging, together with the observation 
that 

m mm 

U -^°,(;) + ^0 = U -^; + ^0 = U ^° + ^0 = Ao, 
we obtain the final expression 

£/(^'t)= -2E E E/H(^)h + ^^-M) 

'^" fc=l 0?^ii;6Ao z=l 

-1 



QtiveKo l<i,k<m 



n 



3.2. Taylor expansion of the energy. We compute in this section the Taylor expansion 
of order 2 of ([14]), viewed as a function on V^^ . To that end, we need to compute 
the gradient and Hessian of Ef at a lattice Aq, respectively at (0,0), and then use the 
approximation 

Ef{H,i) = £^(0,0) + (gradE^(0,0),(H,t)) + ^hess Ey(0,0)[H,t] +o(|| (H,t) f). 

The relevant quantities are given by the following lemma, when / is either an exponential 
fc or an inverse power law ps, which will be the only cases of interest in the sequel. 

Lemma 3.2. Suppose that i^ is such that Aq = [XLi ^'j + ^0 ^s a lattice in R'^. Then 

(1) For an inverse power law psij) = t~^, one has 
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(grad£p,(0,0),(H,t)) = -s £ H[w]\\w\ 



-2s-2 



hessEp,(0,0)[H,t] = s ^ \\zv\\-^^-^ l^l±l(H[w]f -^H^[w]\\wf 

1 9 1 

+^ E 2(s + l)(z(;^(f,--f,^(,-))) -||«;f||f.-i.,(of 

l<i,k<m ) 

(2) for an exponential law fc{r) = e~^^, one has 

(gradE^^(0,0),(H,t)) = -c ^ H[a;]e-^ll«'ll' 






\<i,k<m 

Proof. Using the Taylor expansion of the matrix exponential we write 

exp{H)[w + U - i..,(o] = \\wf + C{H,i) + S (H,t) + o{\\ (H,t) f ) 
where 

and 

5 (H,t) = \\U - t,^ii) f + 2w'H{k - f,,(,-)) + ^h2 [w] . 
Expanding gives 

exp(H)[z(; + f,-f,^(,)]-« = ||a;||-2« A + ^ + ^"j + o(|| (H,t) f) 
in the first case, and 

-cexp(H)[z<;+f,-t (,)] ^ -c\\ivf f , _ ^ ( r ^ ,q\ ^ ^_ tA ^ n( W ( H t\ W^) 



VOJ =e-c|K<'r n-c(C + S) + -C^\ +o(||(H,t) 
in the second one. Then, for a fixed w G Aq, one has to add the terms 

expiH)[w + U - ta.ir)]-', resp. ,-"xp(H)[«;+f,-^(,)]^ 

corresponding to all pairs {i,k). Because aj^ is a permutation, the terms 2w\ti — tg-ii\) 
appearing in C add up to zero, as do the terms 2w'^H{ti — icr^(/)) in S, and the terms 
2w'^{tj — t^j. n\)H[w] that show up in the expansion of C^. Altogether, this leads to the 

formulas of the lemma. n 

10 



There are two noticeable features in the previous calculations, whenever our periodic 
set Aq actually is a lattice : the gradient of the potential energy at Aq, which is a priori a 
linear form in the variable {H, t) G 7id x IR*"^, has a trivial component in the translational 
direction, and its Hessian splits into the sum of a quadratic form in H and a quadratic 
form in t. In other words, when studying local perturbations of energy within the set of 
periodic sets around a lattice, one can consider separatly purely translational moves {i.e. 
with H = 0) and purely lattice moves {i.e. with t = 0). This observation plays a prominent 
role in the results of the next section. 

Remark. The previous lemma extends partly to more general potential functions. For 
instance, one can show, using exactly the same argument as in the proof above, that 
whenever / is a smooth function such that the potential energy E(/, A) is defined and is 
a differentiable function on the space of periodic configurations, one has 

(gradE^(0,0),(H,t))= £ H[w]f{\\w\\). 

4. Main Result 

Using the preliminary computations of the previous section, we can improve the re- 
sults of IICouo6l| . We show that under some rather general conditions, a lattice which is 
locally optimal among lattices regarding energy minimization, is in fact locally optimal 
among all periodic sets. One difficulty in giving a precise meaning to "optimal" or "critical 
point" for the energy, is that a given periodic set admits infinitely many representations 
of type A = [_\^Li U + ^/ foi" various m's and L's. To overcome this problem, we adopt 
the following definition. 



d,m 



m. 



Definition 4.1. Let / be a CM function. 

(1) A periodic set Aq is /-critical if it is a critical point of £(/, A) on V'^^ for every 

(2) A periodic set Aq is locally /-optimal if it locally minimizes E(/, A) on V^^ for 
every m. 

With this terminology, a periodic set is locally universally optimal if it is locally /- 
optimal for any CM function /, or equivalently, due to Bernstein's theorem, for any 
exponential potential fc, c > 0. 

We will allow in some instances {e.g. Theorem [44] below) the wording f -critical for 
a non necessarily CM function /. The least we need is that / is smooth and decays 
sufficiently rapidly so that the potential energy is defined and is a differentiable func- 
tion on the space of periodic configurations. This is the case in particular if / satisfies 
Assumption 2 of Section \T2\ 

Besides the preliminary computations of the previous sections, the main tool we will 
use is the notion of spherical design. 

Definition 4.2. A finite set D of points on the sphere Sr of radius r in R'^ is a ^-design if 

(16) ^T^ f fix)dx = ^ E /(^) 

holds for any polynomial / of degree up to t. 



The following lemma, the proof of which may be found in [Maroi'l or IIC0U06I , gives 



an alternative formulation of the f-design property that will be used throughout this 
section. 

Lemma 4.3 (Venkov IIMaroij Theoreme 3.2.]). Let D be a finite subset of the sphere Sy of 
radius r in R'^ and t an even positive integer. Assume that D is symmetric about 0, i.e. D = —D. 
Then the following properties are equivalent : 

(1) D is a t-design. 

(2) There exists a constant Ct, depending only on r, t and the cardinalty ofD, such that 

xeD 
Our main result may be stated as follows 

Theorem 4.4. 

(1) Let Aq be a lattice, all shells of which are 2-designs. Then, viewed as a periodic set, Aq 
is f -critical for any CM function f, or more generally for any smooth function f such 
that the potential energy E{f,A) is defined and is a differentiable function on the space 
of periodic configurations. 

(2) Let Aq be a lattice, all shells of which are ^-designs. Then, viewed as a periodic set, 

(a) Aq zs locally ps-optimalfor any s > |. 

(b) Aq is locally fc-optimalfor any big enough c > 0. 

Proof. For any fixed positive integer m, we write Aq as an m-periodic set, say Aq = 
\-\T=\ ^f* + ^0- We consider E^ (as in ((T2J)) depending on the particular choice of Lq and t*^ 

to locally study the energy in a neighbourhood of Aq in V^^. We in particular use the 
Taylor expansion of Ef around (0, 0) obtained in the previous section. 

(dJ One has to show that for any CM function /, the gradient of Ef at (0, 0) is orthogo- 
nal to 7id X R'"". Thanks to Bernstein's theorem, it is enough to show it for exponential 
functions fc- For any a > 0, we set 



Ao(a) = < zf G Aq : ||k;||^ = a \ 



These shells of the lattice Aq are assumed to be 2-designs (if non-empty). Using Lemma 
|4.3[ i2J), this is easily seen to be equivalent to the relation 

(17) 2^ ww' = ^ id 

for every positive real number a. In other words, the constant C2 in Lemma |4.3[ [2)) is 

equal to '^^ [j • Observing that H[w] = Tr(Hww^), the expression for the gradient of 
the energy obtained in Lemma [3T2I can be reformulated as 



grad£^^(0,0) = -c £ 



WW e 



t^-c\\iv\ 



2 



OtweAo 
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oi>0 iveAo(oc) 
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-c > > ww'^e ™. 



Consequently, using <^Ty^ we obtain 

d 



(18) (gradE^JO,0),(H,t)) = -c £ ^^^^^^^^^^ Tr(H) = 

a>0 



for any H e Tid = {H e S^ : Tr(H) = O}. 

Thanks to the remark following Lemma [3T2J the previous computation extends readily 
to any smooth function / such that the potential energy E{f,A) is defined and is a 
differentiable function on the space of periodic configurations, since we then have 

gradE^(0,0)= J]/W ^ ww'. 

a>0 weAo{a) 

Thus again (gradE^(0,0), (H,t)) = E«>o ^'^"^'^^"^' Tr(H) = for any H G Tla- 

(tlj To establish local optimality with respect to a given CM function /, it is enough to 
prove that hessEy:(0, 0) is positive definite. By BCouo6l Proposition 1.2], the hypothesis 
that all shells of Aq are 4-designs translates into 

(19) VHGS,(]R), E HH2 = ^l^^((TrH)2 + 2Tr(H2)), 

provided that Ao(a;) is non-empty. By the definition of a spherical design, it is clear that 
a f-design is automatically a f '-design if f > t'. Hence, all non-empty shells of Aq are 
also 2-designs, which implies in particular that 

(20) VHGSrf(R), E H2[H;] = ^^^Tr(H2). 

OtweAoia) 

In case f = ps is an inverse power function, we can plug ( jig) and (J2o|) into the expres- 
sion for hess Ep, (0, 0) obtained in Lemma [3:2} which yields 

hessEp,(0,0) [H,t] = ^|^^^(Lo,s)(TrH2) + -^Y,(t) 

where 

2 



tweAo il<i,k<m 






Unless H is zero, the first term jtjtIt^^^O' s) Tr(H2) is positive because of the assumption 



s > 4. As for Ys(t), we can rewrite it as 



2- 

2 I „,-s-2 



a) l^l</,/c<m 



a>0 09^iy6Ao( 

Since each non-empty shell of Aq is a 2-design, this simplifies to 

ix>Q \ <^ J l<i,k<m 



13 



which is obviously positive for s > j, unless f, = ta^ii) ^^^ every 1 <i,k <m. Given that 
for any pair [i,]) with 1 < z,/ < m there exists k such that 0)^(2) = ; (namely k = crj{i)), 
this last condition implies that f; = tj for all {i,j) and consequently t = mod T. This 
proves assertion (ISaj) . 

U f = fc is an exponential potential, then the same kind of computation as before 
yields 

Tr(H2) ^iu„||2 /^.iu,j|2 /^^ , 1^^ .-c||zH|2 , c 

where 



hess£,,(0,0)(H,tl = ;^ E c||.,|P(c||.,|P-(^ + l)).-** + ^y.(t) 



Oi^iveAo {l<i,k<m ) 

If H ^ 0, the first term ^,\,2) ^Oi^weAo ^ 1 1 ^ I P (^ll^lP ~ (I + -*-) ) e~'^"^" is positive as 
soon as c is strictly greater than j^mTT' where we set 



II 1 1 7 

minAo = min \\w\\ . 

Oi^weAQ 



On the other hand, due to the assumption that all non-empty shells of Aq are 2-designs, 
the expression of Yc(t) reduces to 

a>0 V " / ^<iJc<m 



d 
This quantity is nonnegative for any c > - — -. — — , since it is then a sum of nonneg- 

ative terms, and it is zero if and only if ti = f^- /n for every 1 < i,k < m, that is t = 
mod T. The conclusion follows. n 

As a by-product of the previous proof we obtain that the Hessian of the /c -potential 
energy splits into a sum 

E(^-l)|AoWk-" E ll'i-'..(0f' 

a>0 V " / \<i,k<m 



+ 2 



whenever the 4-design condition is satisfied on each shell of Aq. Here the first term 

pertains to purely lattice changes, and the second term to -purely translational ones. Setting 

c 
y = — , we can rewrite it as 

71 



(21) hess£/^(0,0)[H,t] =y 






14 



TTt/a 



with 

(22) F{y)=i;^(nya-f)\Ao{ix)\e--y^ 

and 

(23) G(i/) = Y, na (nya - (- + Ij j |Ao(a)|e" 

so that, in particular, 

(24) G{y) = -^Hy)- 

With these notations, the assertion that the lattice Aq is locally universally optimal among 
lattices means that G{y) > for all y > 0. If this is the case, equation ([24J) implies that 
f is strictly decreasing on ]0, +00); But then F{y) is positive for any y > since, as we 

already observed, F{y) is positive for any big enough y, e.g. for y > -. — — . In other 

words, we have proved 

Corollary 4.5. A lattice with all of its shells being ^-designs is locally universally optimal among 
all periodic sets, if and only if it is locally universally optimal among lattices. 

This principle applies in particular to D4, Es and the Leech lattice, for which the 4- 
design conditions are well-known to hold (there are basically two proofs, one using 
the properties of the automorphism group, and the other one using theta series with 
spherical coefficients, see IC0U06I for details). Altogether, we obtain 

Theorem 4.6. The root lattices D4, Es and the Leech lattice are locally universally optimal, i.e. 
they locally minimize the f -energy on V^^ for any m and any completely monotonic potential f. 

Proof. Applying the previous corollary, it is enough to check that these three lattices are 
locally universally optimal among lattices. But this is precisely Sarnak's and Strombergsson's 
Proposition 2 in [iSSo6il . The computations on pages 138-139 of their paper show that the 
H-part of the Hessian, 



Tr(H^) ^ ,, ,,2 / „ „2 fd ^ 
H^^ — ^ > c\\w\r { c\\w\r - { - + 1 



e 



OtweAo 
is positive definite for any c > 0. n 
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